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Abstract 



We use dispersion relations to reconstruct, in a model-independent way, 
the formation dynamics of heavy quarkonium from the experimental data on 
e + e~ — * QQ annihilation. We extract a distribution of formation times with a 
mean value for the J/tp, < T j/ip >= 0.44 fm; and for the T, < tj- >= 0.32 fm. 
The corresponding widths of these distributions are given by Atj/^, = 0.31 fm 
and Ar-f = 0.28 fm. This information can be used as an input in modeling of 
heavy quarkonium production on nuclear targets. 



The creation of a heavy quark-antiquark pair occurs at small distances (~ I/tuq) 
and produces compact QQ states which later transform into physical heavy hadrons. 
In the case of quarkonium production on nuclear targets, this evolution can cause 
observable effects Jl]. While these "formation" || effects can in principle be evaluated 
in Glauber-Gribov theory M, in practice this calculation is difficult to perform in a 
model-independent way, since it requires the knowledge of all off-diagonal components 
of the quarkonium-nucleon scattering amplitude. Therefore one often uses a simplistic 
approach, in which the evolution of the quark-antiquark pair is mimicked by a fixed 
"formation time", during which the interactions of the pair are different from the 
interactions of the physical quarkonium; depending on the color state of the pair, the 
interactions of the "unformed" pair can be either suppressed ("color transparency" 
0) or enhanced, if the pair are formed in the color octet state [|. 

In the literature one can find different prescriptions for the formation time r/. A 
still popular viewpoint, for instance, is to assume a universal parameter on the order 
of some characteristic hadronic scale, say, 77 ~ mT 1 . Alternatively, one considers the 
classical expansion of the heavy quark pair and defines formation time as the time 
when the separation of the pair reaches the size of a physical quarkonium state ||. 
Much work has been done also on the quantum-mechanical approach to quarkonium 
formation, where the expansion of a small initial wave packet is controlled by the 
spacings of the bound state mass spectrum ||,0. In this work we address the problem 
of formation time starting from the idea that all essential information about the 
expansion of the wave packet is contained in the correlator of the hard scattering 
operator 0. 

Let us consider the space-time correlator of an operator J, which produces from 
initial state | i > the (QQ) state with certain quantum numbers 

U(x) = (i I T{J(x)J(0)} I i) (1) 

The basic expression which allows the use of experimental data for extracting informa- 
tion from this correlator is the dispersion relation, which in coordinate representation 



10 takes the form: 



1 f 

U(x) = - I ImU(s)D (v^, x 2 ) ds, (2) 

where _ 

D(VSy = ^) = ^-K 1 (V-sr) (3) 



4 7r 2 r 

is the relativistic causal propagator in the coordinate representation; K\ is the Hankel 
function. The expression (§) relates the behavior of the correlator (p]) to experimen- 
tally measurable cross sections for physical processes. For example, in the case of 
(QQ) P arr production in e + e~ annihilation one has simply 

ImU(s) = j^-^ a(e + e- - QQ; s). (4) 
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The physical meaning of eq.(^) is transparent: it represents the correlator as a 
superposition of propagators of physical states, each with the weight proportional to 
the probability of their production in a hard process. 

Thus it is possible to extract information about the space-time evolution of var- 
ious states in a given hadronic channel with fixed quantum numbers directly from 
experimental data. For example, we can define the formation time of the ground 
state with mass m by the time t/ at which the correlator approaches its asymptotic 
behavior 

n(r) ~ r~ 3/2 exp(-mr)] (5) 

note that r = it is Euclidian time. 

To illustrate the notion of formation time in more detail, let us use the following 
simple example - assume that the spectral density (f|) consists of two narrow states 
of identical strength, i.e. 

ImU(s) ~ 5(s-mf) + 5(s-rn%). (6) 
At large time r, the correlator ([!]) will look like 

n(r) ~ r~ 3 / 2 exp(— mir) m^ 2 + m-j/ 2 exp(— (m<i — mi)r) . (7) 
The use of criterion (||) therefore leads to 

1 

r f ~ . (8) 

VTi2 — mi 

Let us now introduce invariant A = t 2 — r 2 and decompose (|3]) in the following way 
& 

D fv^, A) = 6(V\) - 9(\) [Ji(vWA) - iN^VX)] + (9) 

It is clear from (J|) that the formation time 77 extracted from the Euclidian asymp- 
totics (P) of the correlator (fj) will also determine the propagation of the quark- 
antiquark state in Minkowski space, with A > 0. Let us introduce the light-cone 
variables x + = t + z , x~ = t — z, so that A = t 2 — r 2 = x + x~, and conjugate momenta 
p~ = p — p z ~ m 2 /2p z and p + = p + p z ; as can be seen from @, @ and @, in 
Minkowski space the criterion (^) leads to 

l f = Ax + ~ J- ~ 2 2P 2 , (10) 

— m\ 

i.e. the formation length of quarkonium If grows linearly with its momentum p = p z 
in the lab frame. The formation time @ extracted from Euclidian asymptotics (|5[) of 
the correlator is related to the formation length (|Tt]) by the Lorentz transformation: 

'/ — 7 w - T f = — 9 9) (11) 

(m 1 + m 2 )/2 ml — mf 
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where (mi+m^) /2 can be interpreted as a characteristic mass of the wave packet. One 
can recognize this length also as the inverse of the longitudinal momentum transfer 



Ag|| ^ \/p 2 o -ml - \Jpl -m\~ ^ 2 ™ X > (12) 

which accompanies the transition between the two hadronic states at high energy. 
The wave packet can escape from a nucleus of radius Ra before being "formed" if 
If 3> Ra or, equivalently, if 



2 2 

mi — m 



-Ra < 1. (13) 



2p 

This condition was derived many years ago by Gribov ||. 

One can clearly see that the value of formation time is by no means a universal 
constant. It depends both on the properties of the interaction J and on the entire 
spectrum of states in a given hadronic channel. Models using some universal value of 
formation time are therefore misleading. Let us note that eqs.(l-4) indicate that the 
space-time picture of a hard process can be equivalently described in the language of 
the spectrum of hadronic excitations which are formed as a result of this hard process. 
If the operator J is of a short-range nature, then it produces a compact object; Eq.(Q) 
tells us that it is composed of many normal-sized hadronic states. 

Since there exist high quality data on the production of heavy quark states in e + e~ 
annihilation, we will concentrate in the remainder of this Letter on the formation 
dynamics of vector J PC = 1 quarkonium states. In this case the Fourier transform 
of the correlator LT^ = (0|T{ J^(x) J u (0)}\0) has the following familiar structure: 

i J d 4 x e iqx LV(x) = n(q 2 )(q^q v - q 2 g^). (14) 

We use the cross section for e + e~ — ► annihilation, 

Ana 2 

a{e + e~ -> fi + fi~;s) = — — , (15) 

6S 

to express (|]) in terms of the familiar ratio 

„/ \ cr(e + e- -> QQ; s) 

R ( s ) = , + _ + _ r , 16 

resulting in 

ImU(s) = ^R(s). (17) 

Contracting Lorentz indices and using the relation (q^q u — q 2 g^ u )\ f i= u = 3g 2 = — 3<9 2 , 
one can write down the expression (0) as 0, flQfl 



U(x) = -3— / R(s) d 2 D(^,x 2 )ds. (18) 
12ir z Jo 
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Since the scalar propagator satisfies the equation 

at x 2 ^ the final result takes the following form: 

U(x) = — / sR{s)D{^,x 2 )ds. (20) 

AlT z JO 

We parametrize the ratio R{s) in terms of narrow resonances plus a continuum con- 
tribution: 

R{ S )=Y, R i(s)+Rcont{s), (21) 

i 

where 

* W = (2.* + l)» f!'^->g (22) 

is the contribution of the resonance with spin Jj, mass M; and total width IV The 
continuum contribution can be described with a reasonable accuracy as 

Rcont(s) = 3e 2 Q Q(s - s th ), (23) 

where eg is the electric charge of the heavy quark. 

^From (H), (|20|), and (|22"D we calculate the contribution of each narrow resonance 



term 
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where 

Air 



M 2 



(2Jj + -> e+e") (25) 



is the magnitude of the annihilation cross section at the resonance peak yfs = Mj. 
The continuum contribution is 

IWt) = pg-z f°° x 4 K 1 (x)dx, (26) 



57T 4 T" J2M th r 

where the open flavor threshold Sth = AM t h . 

The standard limiting forms of the Bessel function then give the behavior of these 
terms for large and small Euclidean times: 

IT(r - oo) = 3 e-^ (27) 

647r 3 a; 2 v 27TT2 



IWt - oo) = t^L e -^r (2f 



7T2T2 
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„ , . 3cr,T,-iW,- 5 , . 

*< T °> = 64^ < 29 » 

n -(^ - °) = % < 3 °) 

One sees that as r — > oo, the correlator is dominated by the ground state contribution, 
but as r — > 0, the continuum contribution always dominates, independent of the mass 
spectrum. 

We now consider the fraction Fj of a given state present in the QQ correlator 
at a given Euclidian time r, 

F,(r) = ™ (31) 

By definition, Fj(r) < 1, and for the ground state 

Fo(r->oo) = l. (32) 

The formation time ry for the ground state can be defined now by the time at which 
the correlator is dominated by the ground state contribution, or equivalently when 
the relation (|3^) is satisfied with a certain accuracy. 

To illustrate the time evolution of the contents of the QQ correlator, we plot 
in Fig. 1 the functions Fi(r) for J/\l/ and for T. For the charm case, we have 
included the J/ty and resonances, the continuum starting at M t h = Md, and 
also the next three prominent ^(nS) resonances above open charm threshold. For 
the bottom quark case, we have included the T(1S), T(2S), and T(3S) resonances, 
the continuum starting at M t h = Mb, and again three prominent S-wave resonances 
above threshold. To exhibit the effect of the continuum contribution for small r, we 
also show curves for these ratios with the continuum omitted. 

Also shown in Fig. 1 are the discrete values of formation time rj which result 
from the simple estimates in Eq. || using the inverse mass spacing between the ground 
and first excited states. One sees that the rapid increase in the Fi{r) ratios occurs 
generally in the region of values, but of course a range of formation time values 
is involved in the approach towards the ground state dominance of the correlator. 
Operationally, we can then interpret the ratio functions i^(r) as the generalization 
of the form 

Ft\r)^6(r-rf) (33) 

which would be expected if the narrow resonance state had been absent until its 
instantaneous formation at the time ri . We are then led to interpret the derivative 
of the Fi(r) ratios as a continuous distribution Vir) of formation times which occur 
in the evolution of the real correlator, 

Vi{r) = -Ai. (34) 
ar 
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These distributions are shown in Fig. 2, again for the ground states of charmonium 
and bottomonium, and compared with the 5 - function position at the corresponding 
rj°' ) values. One sees that the distributions peak at r somewhat less than these single 
formation time values. The mean values for these distributions are comparable to the 
Tj°\ but the widths of the distributions are also comparable to the mean. Thus in 
evaluating any physical quantity which depends on a resonance formation time, one 
should average this quantity over the distribution in Eq. |33|. 

In principle, information is also contained in the correlator about the time evolu- 
tion and formation of the excited states of quarkonium. However, it is not straight- 
forward how to extract this information, since the procedure used for the ground 
states depends on their dominance at large Euclidean times. We have attempted an 
approximate procedure, which involves omitting in the dispersion integral all states 
with lower mass than the state under consideration. This of course neglects the effects 
of interference in the wave function between the states under consideration and all 
lower states in the wave packet evolution picture, or equivalently neglects all initial 
fluctuations with energy less than the mass of the excited state in the virtual state 
picture. The results for if;', T(2S), and T(3S) are shown in Fig. 3. The general 
shape and parameters are in accord with those for the ground states. It is interesting 
to note that the relatively longer formation time for the ip' is a direct result of its 
closeness to the continuum threshold, whereas one might alternatively attribute this 
to the larger size of the final physical state. Of course, both of these pictures must 
be consistent with the same confinement dynamics, and hence may be expected to be 
related. 

In summary, we have shown how the experimental data on the current correlators 
can be converted to the formation time distributions of the physical states. The 
necessary information in the case of vector heavy quarkonia is provided by the data 
on e + e~ annihilation. The formation time distributions are nonzero in the region of 
time expected from simple arguments involving the inverse bound state spacings, but 
in addition show interesting shapes which persist to large times. Our results could be 
used as an input in phenomenological modeling of quarkonium production in p — A 
and A — B collisions. 
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Figure 1: Formation times for the ground states of quarkonium in e + — e annihilation 
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Distribution of Ground State Formation Times 




Figure 2: Normalized distribution functions for formation times of J/ip and T in 
e + — e~ annihilation 
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Distribution of Quarkonium Formation Times 

Approximate Treatment of Excited States 
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\|/\ <x> = 0.91 fm, Ax = 0.82 fm 

Y(2S), <x> = 0.54 fm, Ax = 0.39 fm 

Y(3S), <x> = 0.59 fm, Ax = 0.42 fm 




Figure 3: Normalized distribution functions from approximate treatment of formation 
times for excited quarkonium states in e + — e~ annihilation 
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